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Abstract 

Recent studies have shown that one-dimensional driven systems can exhibit phase separation 
even if the dynamics is governed by local rules. The ABC model, which comprises three particle 
species that diffuse asymmetrically around a ring, shows anomalous coarsening into a phase sep- 
arated steady state. In the limiting case in which the dynamics is symmetric and the parameter 
q describing the asymmetry tends to one, no phase separation occurs and the steady state of the 
system is disordered. In the present work we consider the weak asymmetry regime q = exp {—[3/N) 
where is the system size and study how the disordered state is approached. In the case of equal 
densities, we find that the system exhibits a second order phase transition at some nonzero (3c- The 
value of j3c = 2tt^/3 and the optimal profiles can be obtained by writing the exact large deviation 
functional. For nonequal densities, we write down mean field equations and analyze some of their 
predictions. 
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I. INTRODUCTION 



Nonequilibrium steady states, wherein the properties of a system are stationary but the 
steady state probabihties are not described by Boltzmann weights with respect to a local 
energy function, may exhibit a number of interesting phenomena absent from equilibrium 
systems: for example, driven diffusive systems Q have generically long-range correlations; 
phase transitions may also occur in one-dimensional non-equilibrium steady states although 
they are of course precluded from one-dimensional equilibrium systems with short-range 
interactions. 

Examples of nonequilibrium phase transitions include the absorbing state phase tran- 



sitions 



2] and boundary induced phase transitions in driven systems wherein a conserved 



current is driven through a finite open system 



a 



Bulk (i.e. not boundary-driven) 



phase transitions may arise in systems with no absorbing states like conserving driven sys- 
tems through the introduction of several species of particles js, 10, 11, 12, 13, 1^ . 



Phase separation has been exhibited in several one dimensional driven systems lla ll^l 
consisting of several species of particles with nearest neighbour exchanges occurring with 
prescribed rates. Models in this class are the ABC model 3|, the AHR model 0,0,0 ^ 
which both contain three species of particles - and the LR model which consists of two 
sublattices with two species each. The phase separation in these models has the striking 
feature of the domains of each species being pure. That is, far away from the domain walls, 
there is zero probability of finding a particle of one species in a domain of a different species. 
This is referred to as strong phase separation. 

An understanding of this phenomenon has emerged through the^act solution of the 
steady state in some special cases of the ABC0 and the LR model|l6|. Even though the 
dynamics is strictly local, in these special cases it has been shown that the steady state obeys 
detailed balance with respect to a long-range energy function. The long-range interaction 
leads to superextensivity of this energy function i.e. the energy of most microscopic con- 
figurations scales quadratically with system size so that the contribution of the entropy 
becomes negligible. Almost all configurations are therefore suppressed and only strongly 
phase separated configurations contribute. Although generally in the ABC and the LR 
models for non-equal numbers of particles detailed balance does not hold and one does not 
have an energy function, the same strong phase separation is observed in simulations 
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In the ABC model there is a parameter q that governs the local dynamics. It describes 
the asymmetry in the rates of nearest neighbour particle exchanges (see section ^J. If 
g < 1 is held fixed then one always has strong phase separation in the thermodynamic limit 
(A^ oo). However for g = 1, in which case the particles exchange symmetrically, the system 
is in a homogeneous disordered state where all configurations are equally probable. In the 
present work we investigate the ABC model in what we shall refer to as the weak asymmetry 
regime. That is, we introduce a system size dependence into q so that an extensive energy is 
recovered (when an energy function indeed exists). It turns out that the appropriate choice 
of q is 

q = exp ( -4 ) • (1) 



Thus (3 is now the control parameter and plays the role of an inverse temperature i.e. 
(3 = 1/T. In this regime an interesting question is as to how the transition from the 
strongly phase separated state to the disordered state occurs as T is varied. Since the 
energy function is long range it is possible to have a phase transition at a well defined 
temperature = (3~^ even though the system is one dimensional (phase transitions in fact 
do occur in one-dimensional systems with algebraically decaying interaction or in the 



Bernasconi model 



2l[ which is an Ising spin model with energy function given by long-range, 



four spin interactions 

In the present paper, we first summarize in section II some known facts about the ABC 
model. In section III, we present our results of Monte Carlo simulations done in the weak 
asymmetry regime for an equal number of particles of each species. These simulations 
indicate that there is a critical value /3c = ^ — 11. In Section IV we derive the exact 
free energy functional for an arbitrary density profile from the expression for the weights in 
the steady state in the case of equal particle numbers. By minimizing this functional we 
obtain the exact value of Pc = ^iiy/S ~ 10.88 . . . and the shape of the density profiles in the 
neighbourhood of the transition. In section V, we investigate the case of arbitrary global 
densities within a mean field approximation. We observe that this approximation turns out 
to be exact in the case of equal numbers of particles. 
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II. THE ABC MODEL 



The ABC model is defined on a Id ring with lattice sites. Each site is occupied by one 
of three types of particles denoted as A, B or C . They exhibit hard-core interaction, i.e. 
only one particle per site is allowed. Neighboring sites on the ring are exchanged according 
to the following rates: 

AB ^ BA (2) 

1 

BC ^ CB (3) 

1 

CA ^ AC (4) 



Thus for q ^ 1 the particles diffuse asymmetrically around the ring and in the case Q' = 1 
they diffuse symmetrically. Note that periodic boundary conditions are implied and the 
dynamics conserves the number of particles. 

This model has been extensively studied in by analytical and numerical means. Let us 
summarise in this section the main results. Starting from an initially random configuration, 
the system coarsens into a strongly phase separated state for any q ^ 1. Thus the steady 
state exhibits long-range order, even though the dynamics is strictly local. We will restrict 
our discussion of g 7^ 1 to the case g < 1 for simplicity. (The case g > 1 can be obtained by 
the transformation g — > 1/g together with the exchange of A and B particles). For g = 1 
particles of all species have equivalent dynamics which results in a steady state in which all 
configurations have equal weight and the steady state of the system is disordered. 

To understand the coarsening dynamics for g < 1, note that domain walls of the type 
BA, CB and AC are unstable to interchanges leading to AB, BC and CA respectively. 
Therefore A particles are driven to the left in a 5 domain and to the right in a C domain 
(particles B and C are also driven in other species domains). Thus the system arrives at 
a metastable configuration of the form A . . . AB . . . BC . . . CA . . . AB . . . BC . . . CA . . . and 
a slow coarsening process, involving the elimination of the smallest domains, ensues. For 
example, the time it takes an A particle to traverse a B domain of length / is of order qK 
Thus the elimination of domains of size / occurs at a rate of order g~' which results in the 
typical domain size growing as Z ~ Int. This growth law which is slower than any power of 



t is referred to as anomalous coarsening 



23|. Ultimately the coarsening process results in a 



strongly phase separated state comprising three pure domains. 
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The general steady state is not, as yet, known. In the special case of equal particle 
numbers Na = Nb = Nc = N/3, however, one can show that the steady state obeys 
detailed balance with respect to a long-range energy function Ti. A configuration of the 
system is specified by the set of indicator variables {Xi} = {Ai, Bi, Ci} which take values 1 
or according to whether site i is occupied by the relevant particle. For example = 1 if 
site i is occupied by an A particle. Clearly, these variables satisfy Ai + Bi + Ci = 1. With 
these indicator variables the energy function may be written as: 

N N-1 
i=l k=l 

and the steady state weights of the system are given by 

P({X,}) = Z^ig^«^-» (6) 



where = ^ denotes the partition function. Note that the energy given in \ll 

differs by a constant from (0). 

In (jnj the interaction between sites i and i + k are both long-range and asymmetric and 
the energy function is superextensive and scales quadratically with system size A^. 

The width of the domain walls in the phase separated state is of order 1/| lng| [15]. So 
for q —>■ 1 the size of the domain walls diverges and the system will be in a homogeneous 
disordered state. Moreover we expect an interesting regime to occur when the width of the 
domain walls is of the order of the domain lengths i.e. 

~ 0{N) (7) 

This yields the weakly asymmetric regime stated in the introduction (0) where the control 
parameter is now T = The steady state weight ()5p6|l also confirms that is the 

natural choice of scaling variable since under this scaling (jHl) becomes 

Pn{{X,}) = Z^' exp {-(3En{{X,})) (8) 

where Zn is a normalization constant (the partition function) and the extensive rescaled 
energy EN{{Xi}) is defined as 

EM{{X,})=n{{X,})/N . (9) 
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III. MONTE CARLO SIMULATIONS 



We have measured by Monte-Carlo simulations the number of nearest neighbour pairs of 
sites occupied by the same species of particles in the steady state. We shall refer to these as 
nearest neighbour (nn) matching pairs. For a completely disordered system (i.e. for /3 = 0) 
the probability of finding a nn matching pair is 1/3 for large A^. As (3 increases, one expects 
this number to increase and to be equal to 1 as /3 ^ oo (i.e. as one reaches the strongly 
separated regime). 

In figure 1, we show the results of our simulations for mjsf defined as 

number of nn matching pairs 1 , , 

rriN = ; : (10) 

system size 6 

1 ^ 1 

= ^ 5Z [ AA^, + B,B,+, + aa+i] - - , (11) 

i=l 

with the occupation variable Xi = Ai,Bi,Ci defined as in Note that because of the 
periodic boundary conditions, site A^ + l is identified with site 1. We see that as increases, 
mjsf seems to be closer and closer to zero as /3 < 11 whereas it seems to have a well-defined 
limit m which depends on /3 for /3 > 11 

m = lim mjsf . (12) 

Note that m o nly contains nearest neighbour correlations as opposed to e.g. the measure 
of order used in [l^ which contains long range terms. Also one could consider the lowest 
Fourier mode of the density profile (see section IV). 




FIG. 1: Parameter m defined by (lllj) for system sizes A = 30 to 240 
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For T < Tc we expect < m < 2/3 with m approaching 2/3 in the hmit T — 0. 
This behaviour can be seen in Fig. ^ For 1/T ^ 11 a crossover from an ordered to a 
disordered state appears which becomes sharper with increasing system size. In Fig. |21 
we plot the specific heat defined as = OEj^/dT. One sees strong finite size effects at 
j3 = 1/T ^ 11. Moreover, the curves suggest that a discontinuity emerges in the infinite 
system hmit which would be consistent with a second-order phase transition. To determine 




FIG. 2: Specific heat for system sizes N = 30, 60, 90, 150 

the critical temperature of the system, we have performed a standard finite size scaling 
analysis [.24j based on the distribution of the parameter rriN defined as in (fTT|). Such an 
analysis has already proved to be effective in the study of non-equilibrium steady states 



25 



26| . We performed Monte Carlo simulations of the model for system sizes = 30 to 
210 over 10^ sweeps in the steady state for each set of parameters. To determine the critical 
temperature Tc we measured the ratio between two moments of the parameter rriN 

Un = l-^rHi ■ (13) 

(In fact, since one does not have to distinguish between positive and negative values of the 
parameter m, the ratio {m?) / (m)^ could as well be used to determine the critical tempera- 
ture.) 

At the critical temperature Tc, Un has a universal value U* . Thus on measuring Un for 
various system sizes as a function of T, U* is the common intersection point of the curves 
and identifies T^.. Our results shown in Fig. |S1 indicate {3^ = l/T^. ~ 10.95. 
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IV. LARGE DEVIATION FUNCTIONAL FOR EQUAL DENSITIES AND THE 
EXACT TRANSITION TEMPERATURE 

In this section we consider only the case of equal numbers of each species 

Na = Nb = Nc = j (14) 
where there exists an energy function given by ()5I9|) 

^ iV-lAf-l 

E^{{X,}) = j;^J2J2^ {B,C,^u + + A,B,^u) (15) 

i=l k=l 

with the occupation variables {Xi} = {Ai, Bi,Ci} defined as in The steady state 

probabilities of the system are given by (jSJ. 

We pass to the continuum limit where pa{x), pb{x), pc{x) are the density profiles of 
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A, B, C particles respectively at position x = i/N. Then the free energy functional 
(or large deviation functional) J-'[pa{x), pb{x), pc{x)] which gives the probability of any 
density profile through 

P[pAix),pBix),pcix)] = exp{-NJ='[pAix), pBix), pcix)]} (16) 
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can be written in terms of the density profiles as 

J^[Pa{,x),pb{x),pc{x)] = K 

+ / dx[pA{x)\n.pA{x) + pB{x)\npB{x) + pc{x)\npc{x)\ 
Jo 

+ (3 dx dz[pB{x)pc{x + z) + pc{x)pa{x + z) + pa{x)pb{x + z)] z (17) 
Jo Jo 

where Pa,Pb,Pc are periodic functions of period 1 and i^" is a normahzation constant such 
that the minimum of JF over all profiles vanishes. The second term on the right hand side 
of (fT7j) represents the entropy of the given profiles while the term proportional to (3 is the 
continuum form of the energy ()15|) . 

It is easy to derive an expression for the optimal profile for Pa{x) by finding the extremum 
of (fT7|) with respect to pa{x) subject to the constraint that 



PA{x)dx = 1/3 . 
One obtains 

Pa{x) = Constant x exp ^~/^ J [Pb{x + z) + pc{x — z)]zdz^ . (18) 
which implies that 

dpA{x) 



-PpAix) (pBix) - pcix)) . (19) 

Similar equations hold for Pb{x) and pc{x) and using pc{x) = 1 — pa{x) — pb{x) one obtains 
the coupled equations 

dpA{x) 



dx 
dpB^x) 



-(3pa{x) {pa{x) + 2pb{x) - 1) (20) 
-(3pb{x) (1 - 2pa{x) - pb{x)) . (21) 



dx 

Clearly one solution of ()20I21|) is Pa{.x) = Pb{x) = 1/3 which corresponds to the disor- 
dered phase, but this extremum is the minimum only in the disordered phase. To test when 
an ordered solution emerges we define the Fourier series of arbitrary profiles as 

oo 

Pa{x) = 1/3 + [a„ exp(i27rnx) + a_„ exp(— z27r?7,a;)] (22) 

n=l 
oo 

Pb{x) = 1/3 + [6„ exp(z27™a;) + 6„„ exp(— z27rnx)] (23) 

n=l 
oo 

Pc{x) = 1/3 + [cn exp(i27rnx) + c_„ exp(— 'j27mx)] (24) 

n=l 



We insert these into ()20l2ip and, anticipating a continuous phase transition, expand to first 
order in a„, &«■ We look for values of j3 for which a solution for nonzero a„, 6„, c„ is present 
and therefore the uniform solution could be unstable. One finds that the uniform solution 
becomes unstable to the nth mode for /3^/3 = (27m)^ so that the first instability occurs at 
pc given by 

Pc = SttVs = 10.882796... . (25) 
Near Pc equations ()2()f21|l can be solved perturbatively in 

P-Pc 



Pc 

One finds to leading order 

Pa{x) = PbI^ + ^I = + ^ 



(26) 



i + (0 2 cos(27r(a; - a;o)) + [£) 2 cos(47r(a; - xo)) + ©(e^/^) (27) 



where Xq can be arbitrary (as an optimal profile remains optimal when it is translated). 
Also, choosing xq = 0, one can show that a„ = a_„ and = for |m| = 0, 1 ... to all 
orders in e and that to leading order in e 

ai~(0''\ «4--(0\ a,r.-Q'^\.. (28) 

The parameter m defined as in 1)121111) becomes in the continuum limit 

/•I 1 

= [Pa + pI + Pl]dx - - . (29) 



and from ()27|) one obtains 



(30) 



Thus the parameter m as defined in ()12lll|) vanishes linearly at the transition. An alterna- 
tive, and probably more standard, choice for the order parameter could be the amplitude Oi 
of the fundamental mode and this would lead to an exponent 1/2. 

We now turn to the calculation of the energy. For arbitrary density profiles the energetic 
contribution to the free energy ()15|) may be written in terms of the Fourier coeffcients as 

~6~^ 2TTin ^ ^ 
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For profiles pb{x) = pa{x — 1/3), pc{x) = pa{,x — 2/3) this energy becomes 

^ oo 

E = - - 3 y sin(27™ /3) . (32) 

n=l 

Thus, using (j^ . near [3c we have ^ 1/6 — 3e/2/5c and the heat capacity —(3'^dE / 3(3 has 
a discontinuity of 3/2 at /^c, consistent with the data of figure |21 

This equal density case is similar to some special cases found in a recent study of the 
dynamical winding of random walks j27| for which the fact that the dynamics satisfies detailed 
balance allows one to write equations for the density profiles of the type and to locate 
the exact transition point. 

V. NON-EQUAL DENSITIES AND MEAN-FIELD THEORY 

We now turn to the case of non-equal densities of particles. A direct consequence of the 
stochastic dynamical rules (jH) is that 

- - - (33) 

and similar equations hold for ^^^^ and We do not know how to solve these exact 

equations, in particular because they require the knowledge of two point functions. One can. 



however, write down mean-field equations 4], by making an approximation which neglects 
correlations (i.e. where one replaces correlation functions such as (Ai^iBi) by {Ai_i){Bi)) 

dt 

-g(A)(5,+i) -g(C,_i)(A) - (5,_i)(A) - (A)(a+i) (34) 
Assuming that the profiles vary slowly with i, we write Pa{x) = (Ai) and 

then keeping leading order terms in 1/A^ and defining r = N'^t in ()36|) yields 

dpA r , , d'^pA 

^ = /?^Mpb-Pc)] + ^ 

t-/3^MPc-p.)].|^ (36) 
dpc nd r , d'^pc 
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One can linearize these equations around constant density profiles 



Pa{x) = rA + Apa{x) , pb{x) = rs + Apb{x) , pdx) = rc + Apdx) 



(37) 



where Apa{x), Aps^x), Apc{x) represent small departures from constant profiles at densities 
va = Na/N, tb = Nb/N, rc = Nc/N (r^, tb, rc are the global densities of the three species 
and of course they satisfy + + = 1). Then one finds that these small departures 
are damped when f3 < given by 

27r 27r 



mf 



rB-rc 



2\l/2 



(38) 



rc = 1/3), the mean field value of 



{2rArB + 2rArc + ^tbTc - r\ 

We see that in the equal density case {ta = tb 
(3c coincides with the exact value ()25|) . Moreover one finds that the solutions of the exact 
equations for the optimal profile (fT^ are steady state solutions of the mean field equations 
(j36|) (as they make the l.h.s. of these equations vanish). So, at least for the equal time 
properties, in this equal density case, (3c and the profiles predicted by the mean field theory 
are exact. 

Near (3ci one can perturbatively find stationary non-moving profiles. The first Fourier 
mode of these profiles is given by 

(39) 

(40) 



Apa{x) ~ V(e) [v^ e2*"(^-^'«) + c.c] 
Apb{x) ~ ^(e 



rc - Ta - tb - ^^Jl - 2{r\ + rl + r^) ^^injx-xo) 



Apc{x) ~ ij{e) 



rB-rA-rc + i\J^ - 2{r\ + r| + 



+ c.c. 



CI ^2i-K(x-xo) _|_ 



(41) 



where e is defined as in (|26|) and 

^(e) 



l-2(ri + r| + r2,) 



a/2 



(42) 



^J2{r\ + r| + rl) - A{r\ + r| + ) 

Analyzing the whole phase diagram predicted by the mean-field equations (j36|) is not an 
easy task. Apart from the constant profile solutions, which become unstable for (3 > 
other (static or moving) solutions might exist in some regions of the phase diagram and a 
full description of the phase diagram would require the knowledge of all the solutions of the 
mean field equations and of their stabilities. 
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Equation ()38p implies that for r\+r'^+r1; > 1/2 there is no second order phase transition. 
However, looking at (jl^ it is clear that the second order transition from the flat profile 
solution to the solution (j39ll42j) should already become first order when r\ + r'^ + r"^ < 
2(r^ + + r^). This is similar to what occurs in a mean-field study of another lattice 
gas[28|. 

At the moment we cannot tell whether the predictions of the mean field theory ()36|) such 
as (jnH|) remain exact in the case of unequal global densities. To try to shed some light on 
this question, we have calculated the large deviation functional to order in a small f3 
expansion. The details given in the appendix show that, in the unequal density case, the 
large deviation functional (|16|17|) becomes 

J^[Pa{x),pb{x),pc{x)] = K+ dx[pA{x)\npA{x) + pB{x)\npB{x) + pc{x)\npc{x)] 

Jo 

+ (3 dx dz z[pb{x)pc{x + z) + pc{x)pa{x + z) + pa{x)pb{x + z)] 
Jo Jo 

- ^(3^ I dx I dz z{l - z) [r^(l - 3rA)pBix)pcix + z) + rB{l - 3rB)pcix)pA{x + z) 
4 Jo Jo 

+rc{l - 3rc)pA{x)pB{x + z)] + 0{l3^) (43) 

We see that for the unequal density case, the large deviation functional is modified and 
terms of order appear which were not present in (fT7j) . This term is the first of a whole 
series in j3. Without knowing these higher order terms, it is neither possible to predict 
the exact value of (3c nor how the equations P9I20I21|) for the most likely profiles would 
be modified. One cannot exclude the possibility that the most likely profile in the ordered 
phase corresponds to moving domains. 



VI. CONCLUSION 



In this work we have investigated a locally driven system of three species on a ring which 
exhibits anomalous coarsening into a strongly phase separated steady state. In the case of 
equal numbers of the particles of each species the steady state obeys detailed balance with 
respect to a long-range superextensive energy function, despite the strictly nearest neighbour 
dynamics. In the weak asymmetry regime where g ^ 1 as in (Q) one recovers an extensive 
energy. In this regime we found a second-order phase transition. 

For the case of equal particle densities we have derived the large deviation (or free energy) 
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functional (I16I17|) for the profiles. As in other examples of non-equilibrium systems studied 
recently [itI ll8|, this functional is non-local allowing phase transitions in one dimension. 
By minimising this free energy functional we have obtained the equations (|19p satisfied by 
the optimal density profiles and analysed the phase transition in the equal density case, in 
particular we found [3c = 27r^/3. 

For the general case of arbitrary particle densities we used a mean-field theory. In the 
case of equal particle densities it turns out that the mean-field solution yields the same 
profiles as the exact free energy optimisation just discussed. An open question remains as 
to what extent mean field theory is valid when the densities are unequal. 

The mechanism for the phase separation that we have studied can be understood through 
the stability of the Fourier modes of the particle densities. In the high temperature phase, the 
system is disordered and the constant density profiles are stable. denotes the temperature 
at which the lowest Fourier mode becomes unstable. As T decreases, more and more modes 
become unstable and the depth of the quench from the disordered high temperature phase 
into the low temperature phase determines the number of unstable modes which can grow 
from the constant profile solution. However, in the steady state one expects only three pure 
domains. How the non-linear evolution of the excited modes combined with the effect of the 
noisy dynamics determine the anomalous coarsening towards the three pure domains is 
another interesting open question. 



VII. APPENDIX: EXPANDING THE STEADY STATE IN POWERS OF THE 
BIAS 

In this appendix, we justify the expression ()4H|1 of the large deviation functional. Let us 
consider a finite system of lattice sites and write the asymmetry q as 

q = e-'^ . 

We are going to show in this appendix that the unnormalized weight P{C) of a configuration 
C in the steady state can be written to order 0^ as 

P(C) = exp[0i?i(C) + 0'i?2(C)] (44) 
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where 



N N-1 



Ri{C) = — — ^ ^ d {BiCi+d + CiAij^d + AiBi+d) 



(45) 



i=l d=l 



and 



R2{C) 



^ ' i=i d=i 

X [iV^(iV - 3NA)B,C,+d + NBiN- 3NB)QA+d + Nc{N - 3Nc)AB,+d] (46) 



In the weak asymmetry regime this leads to the formula ()43|) in the large N limit. 
Let us try to compare two configurations C and C of the form 

C = XiX2...Xi_iAiBij^iXij^2----XN 
C = XiX2...Xi_iBiAi^iXi^2----XN 



(47) 



which differ only by an exchange of a A and a B particle between sites i and i + 1 (the 
notation in ()47|) means that site i is occupied by a A particle, site i + 1 by a -B particle and 
all the other sites are occupied by arbitrary particles). 
If (|45I46|) were true, we would have 



Ri{C') = i?i(C) -3 



Nc 
N 



(48) 



R2iC') = i?2(C) + 



N-l 



Na{N - 3Na) J2 C^+d{N -2d+l) 



N-l 



AN'^{N -I) 

-Nb{N - 3Nb) J2 C^-^d{N -2d+l) + Nc{N - 3Nc) ^ A,+,(iV - 2d + 1) 

d=2 d=2 

N-l 

-Nc{N - 3Nc) B^+d{N - 2d + 1) 

Using the fact that 



d=2 

N-l 



d=2 



(49) 



i+d 



this can be rewritten as 

R2{C') = R2{C) + QaY1 '^^+d{N -2d+l)-QBj2 B^+d{N -2d+l) 



N-l 



N-l 



(50) 



d=2 



d=2 
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where 

^ 3 Nc{N - 3Nc) + Nb{N - 3Nb) - Na{N - 3iV^) 

= 4 iV2(iV-l) (^1) 

- 4 iV2(iV-l) ^^^^ 

and Qc" is similarly defined. Thus exchanging the particles at an AB interface produces 
in ()50|) a difference between a term involving only A particles and a term involving only B 
particles. 

Let us consider now a cluster of m sites occupied by A particles so that C has the form 
C = X1X2 . . . XiAi^i . . . Ai^mXi 

+m+l • • • Xn 

where sites i and i + m + 1 are occupied by particles B or C, so that 

Xi = L' and Xij^rn+i = L" 

where L' is a S or a C particle and so is L". 

There are two configurations C and C" which can be reached from C by single moves at 
the two boundaries of this cluster. A rather simple calculation using shows that 

R2{C') + R2{C") - 2i?2(C) = QA[2mNA - 2 - 2(m - l)N] 

N-l N-1 

+Ql" J2 L\+m+d{N -2d+l)- Qv L'^+diN -2d+l). (53) 

d=2 d=2 

Then summing over all clusters yields 

[MC) - R2iC)] = 2Qa[NI - NNa + {N- 1){Nab + Nac)] (54) 

c 

+2Qb[NI - NNb + {N- 1){Nbc + Nba)] + 2Qc[N^c " NNc + {N - 1){Nca + Ncb)] 

where the sum is over all the configurations C which can be reached from a given C by a 
single exchange and Nab, Nba, Nbc, ■■■ are the numbers of neighboring pairs AB, BA, BC... 
along the chain. 

So far ()48|53j) have been derived assuming that ()45f4fjj) are true. To prove that does 
give the correct steady state weights, one needs to check stationarity, i.e. that 



1 - + y ) [Nab + Nbc + Nca) + A^^^ + Ncb + Nac 



N 



1 - + Y 1 [NBAe^ + NcBe^ + NAce 
+ [^Nabc^^ + Nbcc^^ + NcAC^^) + 0' 5Z [i?2(C') - i?2(C)] (55) 
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This can be checked to order 0^, using simply the fact that for any configuration 

Nab - Nba = Nbc - Ncb = Nca - Nac ■ 

Acknowledgments 

MC would hke to thank the Gottheb Daimler - und Karl Benz - Stiftung, DAAD as well 
as EPSRC for financial support. 



[1] S. Katz, J. L. Lebowitz, H. Spohn, J. Stat. Phys. 34 497 (1984) 

[2] H. Hinrichsen Adv. Physics 49 815 (2000) 

[3] J. Krug, Phys. Rev. Lett. 67, 1882-1885 (1991). 

[4] B. Dcrrida, E. Domany, D. Mukamel, J. Stat. Phys. 69, 667-687 (1992) 

[5] B. Derrida, M. R. Evans, V. Hakim, V. Pasquier, J. Phys. A 26, 1493-1517 (1993). 

[6] G. Schiitz and E. Domany J. Stat. Phys. 72 277 (1993). 

[7] M R Evans, D P Foster, C Godreche, D Mukamel Phys. Rev. Lett. 74 208 (1995) 

[8] P. F. Arndt, T. Heinzel, V. Rittenberg, J. Phys. A 31 L45 (1998) 

[9] N. Rajewsky, T. Sasamoto, E. R. Speer, Physica A 279 123 (2000) 

[10] Y. Kafri, E. Levine, D. Mukamel, J. Torok J. Phys. A35 L459 (2002) 

[11] M. R. Evans, Europhys. Lett. 36 13 (1996) 

[12] B Derrida, S A Janowsky, J L Lebowitz, E R Speer, J. Stat. Phys. 73 813 (1993) 

[13] K Mallick J. Phys. A 29 5375 (1996) 

[14] M. R. Evans, Y. Kafri, E. Levine, D. Mukamel J. Phys. A 35 L433 (2002) 

[15] M. R. Evans, Y. Kafri, H. M. Koduvely D. Mukamel (1998) , Phys. Rev. Lett. 80 425, Phys. 
Rev E 58 2764 

[16] R. Lahiri, S. Ramaswamy, Phys. Rev. Lett. 79, 1150 (1997); R. Lahiri, M. Barma, S. Ra- 

maswamy Phys. Rev. E 61 1648 (2000) 

[17] B. Derrida, J. L. Lebowitz, E. R. Speer Phys. Rev. Lett. 87 150601 (2001) 

[18] B. Derrida, J. L. Lebowitz, E. R. Speer Phys. Rev. Lett. 89 030601 (2002) 

[19] G. Korniss, B. Schmittmann, R. K. P. Zia Europhys. Lett. 45 431 (1999) 

[20] M. E. Fisher, S. Ma, B. G. Nickel, Phys. Rev. Lett. 29, 917 (1972) 

17 



[21] F. Bernasconi, J. Physique 48, 559 (1987) 

[22] E. Marinari, G. Parisi, F. Ritort J. Phys. A 27, 7615 (1994) 

[23] M. R. Evans, J. Phys.:Condens. Matter 14 1397 (2002) 

[24] K. Binder, D. W. Heermann, Monte Carlo Simulation in Statistical Physics, An Introduction, 

3rd ed.. Springer Berlin, Heidelberg (1997) 

[25] O. J. O'Loan and M. R. Evans, J. Phys. A 32, L99 (1999) 

[26] T. Tome and A.Petri J. Phys. A 35, 5379 (2002) 

[27] G. Fayolle and C. Furtlehner, \cond-mat/0211Ul \ 

[28] I. Vilfan, R. K. P. Zia and B. Schmittmann Phys. Rev. Lett. 73, 2071 (1994) 



18 



